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Abstract 

In this paper we study the polarization tensor of photons in a 
QED plasma in the presence of a magnetic field. We do it both 
at vanishing and at finite temperature. We use two different me- 
thods to compute the polarization tensor components IP^ within 
the one-loop approximation. The first starts from the effective 
Lagrangian of the system and relates IT Miy to the thermodynamic 
quantities. The second makes use of the electron propagator in 
an external magnetic field. In this second approach we use an 
imaginary time formalism. These methods give consistent re- 
sults in the first non-vanishing order in the photon 4-momentum. 
Beyond this limit the first method is not applicable and the intro- 
duction of the second can not be avoided. We give some physical 
interpretations of our results. 



PACS numbers: 12.20Ds., 11.10.W, 05.30.Fk. 



1 Introduction 

The investigation of the electromagnetic properties of a relativistic plasma in 
a strong external magnetic field is interesting both from a theoretical point 
of view and for possible applications in cosmology and astrophysics. Large 
magnetic fields are known to be present in neutron stars fl[], where B ~ 10 12 
G, in supernovas 0, where B ~ 10 14 G, and several models foresee very 
large magnetic fields in the early Universe ||. When such strong magnetic 
fields are present the propagation of electromagnetic waves in a plasma is 
considerably modified with respect to the free field case. These strong fields 
can be not treated like a perturbation. Indeed, the tensorial structure of the 
polarization tensor need to be completely reconsidered 

According to the pioneering approach of Fradkin || one of our aims is to 
relate the electromagnetic properties of a QED plasma in a strong magnetic 
field to the thermodynamic properties of the system. This should give the 
reader a physically more transparent interpretation of the results that par- 
tially are present already in the literature, as well as of the new results that 
we obtain here. 

All the thermodynamic properties of a QED plasma can be derived start- 
ing from the effective Lagrangian of the system. Assuming B(x) = (0,0, B) 
the Lagrangian can be written 



eff 



C eS (B)+C eS (B,T,ft), 



(1) 



where the vacuum contribution at the one loop level is |j] 




whereas the one-loop matter contribution is 



£ eS (B,v,T) 



1 \e\B 
/3 2^" 



' pv 

/ dk 3 {ln[l + e 




Here 




(4) 
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is the energy of the n-th Landau level at tree level. The sum J2 n =6 = \ Sn(2— 
S n o) takes into account the double spin degeneration of the Landau levels with 
the exception of the lowest. Equation (|3]) follows from a simple phase space 
consideration []. In fact, due to the coalescing of the transverse momentum 
states into those of a two-dimensional harmonic oscillator, the number of 
available states for any given value of n is 

To C ^ ( 2 ~ *W dka . (5) 

One of the most important effects of the modification of the electron 
phase-space is the "pair generation" that takes place when magnetic fields 
larger than B c = m 2 /e are applied to the plasma. We use quotation marks in 
order not to confuse the reader about the meaning of pair generation in this 
context. Clearly, a constant magnetic field can not induce pair generation 
from the vacuum. However, things can be different at finite temperature 
and/or density. Indeed, the number density of electrons+positrons is ]7] 

\e\B ~ f 1 1 I 

/ rffc 3 1 1 i „3(Er,(k*)-u) + 1 I „3(Er,(k*) + U,) I ' W 



n 



2tt 2 ^ 7-oo I 1 + e P{En{k 3 )-li) 1 + e 3{E n {k 3 )+v) 

n=0 

It increases roughly linearly with B when B ^> B c and T <{\e\B) 1 ^ 2 . This 
phenomenon can be understood as a shift in the equilibrium of photons and 
pairs ||: the equilibrium of the process e + e~ <-> 7 moves to the left owing to 
the growing of the number of available states that electrons and positrons can 
occupy in the lowest Landau level. Thanks to the amplification of its phase- 
space this level is practically the only occupied level when eB ^> T 2 . Since 
B modifies the density of the charge carriers, it is reasonable to expect that 
the electromagnetic properties of the plasma have to be affected by strong 
magnetic fields. The charge neutrality of the plasma is preserved, provided 
that fi = 0, since the increasing of the electron and positron energy densities 
balance each other. Indeed, the charge density is 

dC cS (B,T^) 

M = ^ = (7) 

\e\B ~ roo t r 1 1 



2tt 2 

n=0 



^ r ( 1 1 1 

^ J-00 3 1 1 + e^ E ^~^ 1 + e /3(^«(fe3)+M) J 



1 Although Eq.(||) can be obtained from the fermion propagator as done in refs. 
we prefer to adopt here this simpler physical interpretation. 
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The magnetization plays an essential role in determining the screening 
properties of the plasma. Like the energy density, the magnetization M can 
be obtained starting from the effective Lagrangian 

dC dC° s (B) , djC°*(B,T,n) 
M -dB~^B~ + dB • (8) 

For example, the matter contribution to the magnetization in the limit T = 
is (we refer the reader to ref. || for a more general expression) 



M(B,0,/i) 



e / r~ n : — ~ ~ , . „ u + vw - m'- x 



/xy/i 2 — m 2 — 2eBn — (m 2 + AeBn) log 



i//i 2 — " nn ' 1 



27r 2 a V y/m 2 + x J 

(9) 

Note that in contrast to what happens for p, a component of the magne- 
tization depends on the vacuum part of the effective Lagrangian. In other 
words, the vacuum polarizes in strong magnetic fields. The relative strength 
of the vacuum and matter polarization has been evaluated in ref. ||. The 
vacuum contribution to the magnetization exceeds the matter contribution 
if eB 3> T 2 , \fi 2 — m 2 \. We do not consider here the other vacuum contribu- 
tions to the polarization tensor since they are already widely reported in the 
literature (see e.g. ref. fllPf). 



In the next section we will show that the components of the polarization 
tensor can be related to the thermodynamic quantities p, M and the magne- 
tic susceptibility x — OM/dB. In section 3 we write the general expression 
of the polarization tensor in terms of the fermion propagator in an external 
magnetic field. In section 4 we apply this expression to get the components 
of IP 1 ' in the static limit and verify that they coincide with the results ob- 
tained in section 2. In section 5 we give examples of computations of some 
components of beyond the static limit. Finally, section 6 contains our 
conclusions. 
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2 The screening operator in terms of £ eff 

The polarization operator for a QED plasma is defined by 

IV x, x) = i — . = i -r— TT r- * / .\ . io 
M v y <5A„(x) 5A,(x) &A M (a/) 

where T cfr = / d 4 x £ eS is the effective action. If the plasma has a non- 
vanishing chemical potential the tree level Lagrangian is 

£ = -lF^F^ + ^-e4- 7o //-m)V> • (11) 

The form of Eq. (|10"D remains unchanged if we replace A{x) with A(x) = 
(A + i/x, A). Note that only this combination of the chemical potential and 
the vector potential is physically meaningful [[11]]. In fact, whenever A is 
changed, the polarization of the plasma rearranges so to fulfil /i+eA = const, 
at the equilibrium []E|. Thus, /x is determined by the charge distribution of 
the plasma, and its value may depend on the position. It is then convenient 
to work in the gauge Aq = in order to have a well defined A(x) = fi = const. 
for the whole plasma. 

As a first application of the above formula we compute the limit for x — > x' 
of II 00 . We then have 

n (B,T^) - e — = e^= (12) 

\e\ 3 B(3 ~ f°° ( 1 1 



£ / dk 3 

* Jo 



2 ( P{E n (k z )-ji) 



(2vr) 2 ^J ^cosh 2 ( ^" ( 2 fc3)+ ^ ) cosh\ , 

where, assuming a uniform B and working in the Euclidean space, we have 
used rf at = -i(3V£(B, T,/i). 

It is worthwhile to compare this result with the corresponding quantity 



obtained for B = [13 



Keeping in mind Eq. (J5|) it is evident at a glance that the effect of the magnetic 
field on II 00 is completely mediated by the modification of the electron phase 
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space. The increase of the screening power of static electric fields when 
overcritical magnetic fields are applied to the plasma can also be interpreted 
as due to the increase of the number of pairs (see Eq. @) as we anticipated 
in the introduction. In this sense the magnetic field plays a role analogous 
to that played by the temperature. 

Let us now consider the II * components (i = 1,2,3). In the same gauge 
as we used above, we have 



n 



Oi (•£) % 



5 5 



^ r 



eil' 



ic e 



ijkdft 



dp 



5 A 5 A 1 '""-"SB^x') 
where we used the chain rule for the functional derivative to get 

identically and n i 



(14) 



5 



5A i 

if the momentum 



Whereas H03 = identically and n i = II02 

insertion is zero or parallel to the field, the transversal components receive 
a non- vanishing contribution if the momentum insertion is perpendicular to 
the field. More precisely, the Fourier transform of (14) in terms of the photon 
momentum components pi is 



n 



0/ 



dp o 



1,2 



(15) 



while rijo = n^. These terms are related to the to the interaction of the 
net charge of the plasma with the magnetic field and contribute to the Hall 
conductivity. For the spatial components we get the general formula 



n 



ij 



kl ran 



PkPr, 



B n BA 1 d , B n B { d 2 



B 2 J BOB B 2 OB 2 



C 



cff 



This gives 



n 33 


2 i 

= P2 B 


M ■ 


n n 


2 1 

= P *B 


M + p 2 2 x ; 


n 22 


2 1 

= 


M + P lx ; 


n i2 


= n 2 i = 


= P1P2 x ; 


n 13 


= n 3 i = 


= -PaPi-g M 


n 23 


= n 32 : 


1 

= pm^ M 



(16) 

(17) 

(18) 

(19) 
(20) 
(21) 

(22) 
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It is worthwhile to note the symmetry of the spatial part of the polarization 
tensor. This agrees with the tensorial structure obtained in ref.0. 



3 The electron propagator in an external ma- 
gnetic field 

Although the results of the previous section have been obtained without any 
use of the fermion propagator S, a treatment addressed to determine the 
components of the polarization tensor beyond the leading order in momenta 
can not leave S out of consideration. The fermion propagator in a constant 
magnetic field is given by [[L4]] 



S(x, x') 



E 

n=0 



(27T)' 



where 



uj 2 — k\ — m? — 2eBn + ie 
S(n;u,k 2 ,k 3 ) = 



S(n;u,k 2 ,k 3 ) 



(23) 





-(u- k 3 )I n 





i\ / 2eBnL n _ 1 



i\f2eBnI n _ 1<n 



\ -iy/2eBnI n -i, r . 
and 







(u + fc 3 ) I n -l,n-l 

In,l — In;k 2 ( x )Il;k2\ x ) 



where 



■M x ) = — 



eB\ l ' A 



7T 



1 



2 n/2 



-iV2~eBnI njn -i \ 

-(iO — ^3)/ n -l,n-l 



m/ n _ 1(n _i / 
(24) 

(25) 



GB[x-£ l (26) 



H n are the Hermite polynomials and we have chosen the 7-matrices in the 
chiral representation. The propagator is obtained by solving Dirac's equation 
in a magnetic field. Some useful identities are 



In;k 2 (x)In;k 2 (x f ) = 6(x - x') 



(27) 



n=0 
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and _ 

/oo 
dxI n . M (x)Ii- k2 (x) = 8 n>l . (28) 
-oo 

In terms of S(x, x') the polarization tensor is 

W v {x) = e 2 J dt'dx'dy'dz'Tr (^S(x, x')YS{x', x)). (29) 

Before we pass to the computation of the components of IT Miy beyond 
the static limit we will show how the results of the previous section can be 
reproduced using Eq. fl2"9"|) . 



4 The static limit 

Whenever the momentum insertion from the electromagnetic field into the 
fermion loop is vanishing the polarization tensor can be written 

x " dudk 2 dk 3 1 au 



n^= E J 



n,l=0 ' 



(2tt) 3 (uj 2 -kl-m 2 - 2eBn) (u 2 -k\-m 2 - 2eBl) 

(30) 

where 

tt^ = J dx' Tr ( 7 ^(n; u, fc 2 , k 3 )j v S(l; u, fc 2 , fc 3 )) (31) 
It can be calculated, for example, that 

vr 00 = 2{u 2 + kl + 2eSn + m 2 ) (J n>n + / n _i >n _i) <f,, n (32) 

vr 33 = 2{uo 2 + kl - 2e5n - m 2 ) (J n>n + / n _i, n _i) (33) 

vr 11 = 2(u; 2 - fcf - m 2 )(/ njn ^ in+ i + 4-l,n-l^,n-l)+ 



2e5 Jn(n + l)J n _ l n+1 <y, n+1 + Jn(n- l)J n ,n-2*i,n-i (34) 



7T 11 differs from 7r 22 only in that the last term has the opposite sign. As 
expected II 00 and IT 33 correspond to zero angular moment transfer between 
the photon and the plasma, implying the selection rule An = 0. Instead, 
transversally polarized waves can induce transitions between different Landau 
levels. 
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4.1 The T = finite \i contribution 

Let us consider the T = contribution for the various components. We have 

^m,^ s 8|e| 3 5 ^ r duo f°° „ cj 2 + k\ + 2e5n + m 2 , . 
V W (2vr) 2 ^Jc2mh (u 2 - kj - 2eBn - m 2 ) 2 V ; 



To obtain this, perform the &2 integration using (28)- It follows that 

8\e\ 3 B ™ r du , 7 ( 2u 2 1 

22f7r-l dk 3 



(2tt) 2 ^ J 2ni Jo {(ujz-k 2 - 2eBn - m 2 ) 2 oj 2 - k 2 - 2eBn - m 2 < 



8|e| 3 5 ~ / d /-v 7 ^^ „ /-= fy /tf-m*-x dk 



Eh-/ ^ 3 \Mf + ™ 2 + ^ 



(2tt) ^ ^dx-A) Jo 2 v /A;l + m 2 + x 

7T 2 ^ V^-m^- 26715' 1 ' 

71=0 

where x = 2eBn. As can be seen from above, the only no n- vanishing contri- 
bution is coming from the point where the uj = E double pole is crossing the 
integration contour. More about the contour C and the methods we use to 
perform the integrals can be found in the Appendix. 

There are two regions where the expression simplifies. One is the limit 
B — where we get 

2 /X\//i 2 — ™ 2 



e (JVM — 771 



n 00 (0,0,/i) = ^ v ^ 2 - (37) 
the other is 2eB > fi 2 — m 2 where 



um ^ J ^7F^ ' (38) 

Eq. ( |38| ) can be understood in terms of the modified relation between the 
charge density and the chemical potential of the plasma. In fact, in this limit 
Eq. (H) becomes 

eB 



PO) ~ V -™ 2 (39) 
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and Eq.fl38|) can easily be reproduced using II 00 = e 2 dp/dfi. Eq. fl3T|) di- 
sagrees with the result of refs.^]- In fact, the authors obtain zero for 
n 00 (0,0,/i). On the other hand, our Eq.(|5()|) agrees completely with the 
result of ref . |T5| . 

Let us now consider the spatial components and verify that they are zero. 
From (1331) we have 



^•H/r. n \ S\e\ 3 B ™ r du f°° „ u 2 + k\ - 2eBn - m 2 . . 

n 33 (5,0,/i) = J-l £ f — dk 3 J -3 = (40) 

(2n) n=d J 2m Jo [u 2 - k% - 2eBn - m 2 ) 

Again it is important to take the contour crossing into account in order to 
obtain the correct result. 

I1 11 (S, 0, n) and LT 22 (B, 0, /i) are simply computed. The last term of (134]) 
is killed by the k<i integration and we are left with a vanishing result. Hence 
we conclude that 

u n (B, o, n) = n 22 (5, o, fi) = n 33 (s, o, fi) = o (41) 

Thus static magnetic fields are not screened in the present limit. These 
equations clearly remain valid when B — > 0. Again we disagree with the 
conclusion of ref. [1^] where a magnetic screening length proportional to the 



Fermi momenta was obtained at B cxt = 0. The disagreement that we ob- 
tain here and for n 00 (0,0,/i) suggests to us that the contributions from the 
crossing of the poles through the contour C was ignored by the authors of 



ref. P| 



Proceeding similarly to what we have done for the other components we 
easily deduce that IT - 7 (5, 0, /x) = if there is no momentum insertion in the 
fermion loop. 

4.2 The finite T contribution 

At finite temperature we must add two more contour integrals: 

v ' n=0 N 



+ 



-ioo+fi+e 



~ai rllLO, k,n) + q> l(oj,k,n)\ 

2m l + e^-v) v ' ' ; J c 2m v ' ' ') 
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where 



I (id, k, n) 



uo 2 + k 2 + 2eBn + m 2 
(uj 2 — k 2 — 2eBn — m 2 



(43) 



The integrals are performed according to the prescription given in the 
Appendix with the result 



n 00 (B,T,v) 



^(3 



47T 2 



E 

n=6 



2 (§jE±jx)\ 



cosh 



+ 



1 



cosh 2 fi&A 



in agreement with (IE 



ref. 13 in the limit B 



(44) 

. Although here the result coincides with that of 
0, this can easily be understood since the contribu- 
tions from the pole crossing through the contours C and C± cancel in this 
case. Similarly one can verify that the spatial components remain zero even 
at non-zero temperature. The behavior of n 00 as a function of B is shown 
in Fig. 1. As we anticipated in the introduction, we can see from this figure 
that the screening properties of the plasma are sensitive to the magnetic field 
only when eB 3> T 2 . Thus the effect of even a large magnetic field can be 
neglected if this condition is not fulfilled. 



o 
o 




Fig.l. II for three different temperatures. From top, T = 10m e , T = m e and 
T = 0.1m e . 
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5 Beyond the static limit 



The application of an external magnetic field breaks the isotropy of the 
plasma. It is thus understandable that the electromagnetic response of the 
plasma will depend on the direction of the the momentum insertion with 
respect to the field orientation. In this section we will give some sample 
calculations when the photon has non-zero momentum or energy. 



5.1 p || B 

We start by considering a momentum insertion parallel to the magnetic field. 
Put k' 3 = k 3 + p 3 where ^3 is the momentum insertion. IP 1 ' is then obtained 
by replacing k 2 with k' 2 in one of the propagators and replacing k 2 with k 3 k' 3 
in 71^ in Eqs. (|32| - |34"|) . To compute II 00 we use 

ti r ^ uj 2 + k 3 k 3 + 2enB + m 2 

h - n) = W - W) ^ - eMF) (45) 

as the integrand in Eq.(35). Performing the uj integration over the contour 
C (see the appendix) this becomes 



m 2 + 2eBn j fy/v 2 -m 2 -2eBn dk 3 . 

(P3 -> -p 3 ) 



2p 3 \J-y/^-m?-2eBn (2k + p 3 )\/k 2 + m? + 2eBn 
Expanding to second order in p\ we find 

n 00 ^, o, fi\ P3 ) = u 00 (b, o, flips = o)+ 



(46) 



2eB 



P 3 \ e \ B I » v 



vr 2 ^ \12(/i 2 -m 2 -2eBn) 3/2 6y ' fi 1 - m l - 2eBn(m 2 + 2eBn) 

+ 0(pl) (47) 

where U 00 (B,0, fi) is given by Eq.(36). 

At T ^ we have to manage two simple poles in the complex cu-plane. As 
discussed in the Appendix, the integration over C + splits in two integrals each 
of them corresponding to the two possible positions of the positive-energy 
pole: E > fi or E < fi. A third integral, over C_, receives a contribution only 
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from the negative-energy pole, always contained in C_. An T-independent 
term coming from the integration over C + and C_ cancel exactly with the 
result of the integral over C given by Eq.(^Bj). Thus we get 

n (B,T^\p 3 )-—^^ x 

f°° dk ( 1 1 ^ (4&) 

J-oo (2k + p 3 )Vk 2 + 2eBn + ^ VI + e^ w ~^ + 1 + e ^ w+ ^ ) ' ' N ' 

Previously we found that U 33 (B,T, /x) = at zero momentum. At non-zero 
p 3 we have to integrate 

uj 2 + k 3 k' y — 2enB — m 2 

(49) 



(u 2 — E n (k 3 ) 2 ) (uj 2 — E n (k' z ) 2 ) 



h 

over k 3 and uj. The residue at uj — E n is . Clearly, II 33 = even at 

2E n p 3 

non-zero p 3 after integrating over both signs of k 3 . This agrees with Eq. (p~7|) . 

It is also worthwhile to check II 11 with non-zero p 3 . Using Eq.(|34]) we 
find that we have now to manage the integrand 



uj 2 — k 3 k' 3 — m 2 



(uj 2 -kl-m 2 - 2neB) {uj 2 - k'£ - m 2 - 2(n + l)eBl) ' 
After the integrations we get 



(50) 



n 11 ( J B,o, / ib 3 ) = § 



\e\ 
27 2 



n=0 



( fiyj /i 2 — m 2 — 2eBn 



[m + AeBn) log ^=^= 

ym? + x 



0(pi) (51) 



that agree with Eq.fllSD as the reader can check keeping in mind Eq.@. The 
same result as (|5T| ) is found for II 22 (B, 0, fi\p 3 ). 

Let us now move on to IT - 7 . Using (|3l|) we find that 



7T 



03 



2uj (k 3 + k' 3 ) (I n>n + 

In— 1,11— 1 

)S n>l (52) 
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This clearly gives zero even for p^ ^ 0. In fact, the integrand can be rewritten 

as 

oj / 1 1 \ 

P3 Vtu 2 — k 2 — 2eBn — m? to 2 — A;' 2 — 2eBn — m? J 
Shifting the k' in the second term gives zero after integration. We leave to 
verify n i(-B, 0, /j.|p 3 ) = U 02 (B, 0, n\pz) = as an easy exercise for the reader. 

5.2 p±B 

The case with momentum transverse to the magnetic field is slightly more 
complicated. Of special interest are the II 01 (i = 1,2) components at non- 
zero pi which give through 

dll 0i (v) 

o» = i^r^ 1 (54) 
the plasma conductivity. Using ([y]) we find that 

■•■/wvan,,^^ (55) 

At non-zero P2 we have to calculate 

dx'dy'dk 2 dk' 2 e^^-^)- 4 ^^^')^^-^)^^^^)^^^^)^^^^)//^^^) = 



eB J dxdx' I n;0 (x)I n ,. (x')Ii m (x)Ii> ;P2 (x r ) (56) 
This can be found to be 

for n > I and n' > /', where the L's are Laguerre-polynomials. Expanding to 
order p 2 we find, in our case, 

- V2^Bp 2 (I5 nj+ i + I8 n j) + (n -> n + 1) (58) 

and to order pi 

[3l 2 + l) 5 nHl - (l 2 + l) 5 n , l+2 + (3/ 2 - l) 5 n>t - (/ 2 - l) 5 n ^) 



Pa 



2V2eB 



13 



+ (n-*n + l). (59) 

For larger powers of p 2; larger jumps between the Landau levels are allowed. 
We get 



2 2 
— m 

2eB 



U°\B, 0, An) = ' £ ' f ^ -m»- 2e5n e5n 



7r 2 . \ v v 'u 2 — m 2 — 2eBn , 



[' 



Y. UJ? - m> - 2egn - - \ + 0(rf) . 



(60) 

n 2 is given by the same expression with p 2 replaced by p\ . Eq. (pDf) confirms 



the result for the conductivity of |Lo] and extends it to second order in mo- 
mentum. Note that the higher order correction vanishes for 2eB > fi 2 — m 2 . 
The expression for the conductivity given in ref.[16 , on the other hand, do 
not agree with our result or the result of fl5 l. 

We are also able to extend our result to finite temperature. 



n=0 



j pr f 00 / A \ poo / 1 1 



4n 2 B^\ dBjJo 6 yi + eW-ri 1 + e^ E+ ^ 

+ 0(pl) (61) 



The first order term can again be checked using (15). 

We may also consider II 00 and II" at nonzero p\ and p 2 ■ The calculation 
uses the same method as above. After expanding the Laguerre-polynomials 
and performing the contour-integrals, we can verify Eqs. (17-19). On the 
other hand we find 

n 00 (5,o,/x|p 2 ) = n 00 (5,o,/i,p 2 = o)+ 



e 2 p\ ~ / r- — — - 2 n + vV 2 -m 2 -x AeBnfi 

2 t /ly fj, — m z — 2eBn + m log ■ 



2v 2 eB^ n V ' Vm 2 + 2eBn Vv 2 - m 2 - 2eBn 

+ 0(pi) (62) 

Together with ( |47|) we then have the complete expression to second order in 
momentum. In this case we can not check the result using section 1. 
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5.3 Energy insertion 



Let us now consider the case in which the photon energy is different from 
zero. In this case the energy of one of the fermions in the polarization loop 
will be shifted to u' = u> + p where the photon energy p is assumed to be 
real. To compute IT 00 (.B, 0, /x), the integrand to replace the one in Eq.(35) is 



I(u,k 3 ,n) 



ujuj 1 + k% + 2enB + m 2 



E n {k 



: 3 J 2 ) (a/ 2 - E n {h 



(63) 



Performing the (^-integration we find 

n 00 (B,o,^\p ) 



/r..2 



\e\ 3 B 

p lT 2 



(m-po) 2 -™ 2 



V 



y /i 2 — 2enB — m 2 

n=6 



2\e\ 3 B 



2 2 
fi — m 

2eS 



2eS 

£ 

n=6 



2enB + 



\ 



(/i — po) 2 ~~ 2enB — rri 



J 



(5,0, ^+—5- ^ ^°2(^ 2 - 2en5 - m 2 ) 3 / 2 +Po (/x 2 - 2enB - m?f'\ 

+ O(p 3 ) . (64) 

At T 7^ we have instead 

8|e| 3 5 ^ r r°° / 1 1 



n 00 (s,r )A # ) 



(2vr)% 5 



dk-> 



3 ^ e ^(E-M) + 1 e /3(S+M) + i 



fife 



1 



1 



z j3(E-(p,+p )) _)_ ]_ e /3(S+(M+Po)) _|_ 1 



(65) 



6 Conclusions 

In this paper we have studied the polarization tensor of a QED plasma placed 
in a magnetic field. Although this has been the subject of several studies, 
some discrepancies, even in the simplest limit B = ||13|| , called for a more 
careful analysis. Furthermore, recent results concerning the thermodynamics 
of a QED plasma in a strong magnetic field allow a new, physically more 
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transparent, interpretation of old results as well as of the new results that 
we obtain here. 

We have used two different methods. The first relate the polarization 
tensor directly to the thermodynamic quantities of the system. The second 
makes use of the fermion propagator in an external magnetic field. These 
methods give the same results in the first non- vanishing order in the photon 
4-momentum. Beyond this limit the first method is not applicable and the 
introduction of the second can not be avoided. For instance, this is needed 
if one wants to find the dispersion relations for the electromagnetic waves. 

Some of our results have a simple physical interpretation. Static electric 
fields are screened by the plasma in strong magnetic fields (eB > T 2 , fi 2 — m 2 ) 
more effectively than in the free field case. This can be understood since static 
electric fields are screened by the charge rearrangement in the plasma. At 
T = no thermal pair production is active and the screening can be achieved 
only in the presence of a charge asymmetry. In strong magnetic fields the 
relation with the charge asymmetry is modified. Indeed, the charge density 
becomes proportional to the magnetic induction when eB is larger than the 
Fermi momentum squared, /i 2 — m 2 (see Eq. (|39|) ). Consequently, IT 00 grows 
linearly with B (see Eq.(|38|). 

At T ^ 0, as we discussed in the introduction, the thermal pair production 
is amplified by large magnetic fields. Then more charge carriers are available 
to screen static electric fields. The reader can see in Fig. 1 that II 00 starts to 
grow linearly in B, once all electrons and positrons have dropped in the lowest 
Landau level. In this static limit the anisotropy induced by the external 
magnetic field can not play any role and the definition of an electric screening 
length A~ 2 = m 2 el = II 00 (p = 0, p — > 0) |TJ] is still meaningful. 

Static magnetic fields are not screened at all. In fact, macroscopic spatial 
currents can not be obtained if p = (see Eq.|16|). This conclusion can not be 
modified in presence of strong external magnetic fields neither at vanishing 
nor at finite temperature. However, magnetic screening is achieved if p ^ 
and Eqs . ([T7| - |22D show how this is related to the magnetization and magnetic 
susceptibility of the plasma in a non trivial way. 

We also computed the electric conductivity of the plasma. Our result 



confirms that of ref . |15 1 , but we improve the calculation by including second 



order terms in the momentum expansion. On the other hand, we disagree 



with the result of ref. 16 



The determination of the complete dispersion relations for electromagne- 
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tic waves propagating trough plasmas in strong magnetic fields is beyond the 
purposes of the present paper. One would need to take into account the full 
tensorial structure of the polarization operator as done in Refs.[||. Using 
our results according to the prescription of Refs.jp, would provide physically 
more accessible information about the propagation of electromagnetic waves 
trough a plasma in presence of strong magnetic fields. 
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Appendix 

Let us give some details on how the contours are specified. At finite temper- 
ature the integral J du is replaced by a sum over u = ^f(2n + 1). This can be 
accomplished by integration along a counter clockwise contour around the 
poles of tanh e.g. 

1 f 7T Z 

£ f(2n + l) = -fdz tanh —f(z) (66) 
Deforming the contours gives 

Loo 2^ /M L W e 2^ /( "W) + l + (6?) 
rioo+n-e duj 1 f dUJ 

/ t; ■f\ U} ) - 0(ul - u ) i i + f t^—J^) 

J-ioo+v-E 2m e W + 1 Jc 2m 

The first term is the divergent vacuum contribution, the second and third 
terms are non-zero only at finite temperature and the last term is the T = 
non-zero \i contribution. The contours are depicted in fig. 2. 
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r" 



A A 



c. 



Fig. 2. Contours for T = on the left and for T ^ on the right. 



It is worthwhile also to give some details about how we performed the integrals. 
Mainly we wish to direct the attention of the reader to the crossing of poles through 
the contours of integration. 

We start by considering the T = finite-density contribution. This term has 
the form 



duo 
C 2ni 



dkf(uj, k) =Y^ Hm^-^-B. 



dw r 



dkf(w,k)(w-Ei)?0(jM-w) 



(68) 

where we have assumed that the function f(u>, k) has poles u = Ei on the real axis 
of order mj. The contour of integration C is depicted in Fig. 2. The result of the 
integration is zero if the poles lie outside the region delimited by the contour C. 
This justifies the ^-function on the left side of Eq.(|68[). The effect of the ^-function 
is simply to modify the effective integration limits in the case the pole is first-order. 

More interesting, when the pole is second-order (we do not come through 
higher order poles in our one-loop computations) a new term appears containing a 
delta-function coming from the derivative of the 9. This term corresponds to the 
crossing of the pole trough the contour of integration. The 5 kills the ^-integral 
and leaves us with a term of the form like that in the result of Eq.(36). 
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The finite temperature contribution has the general form 

/ dk f{uo,k)—, , + (69) 

/ dk f(u,k) — ^7 ; h <p / dk f(u!,k) 

When the poles are first order the first two, temperature dependent, integrals 
split in three non- vanishing parts. The first part corresponds to the integral over 
the contour C + (see Fig. 2.) for values of k such that the pole u = E is inside C+ 
(E > /j,). It can be written 



1 P \^ fJi ~ — X — VYlP' 1 

ik l^jm^- Ljs^ ik ,{E ' k) im-WTi ; (70) 

the second part corresponds to the integral over C_ when only the pole u = —E 
is inside C_ 

dkf(-E,k) 1 . (71) 

and the third one comes from the contribution to the same integral when k is such 
that the pole to = E is inside C_ 



\J fj? —x—m? 1 



It is then straightforward to verify that the sum of the second integral in ( |70|) and 
the integral in ( [72]) i s temperature independent and cancels with the integral over 
C in ([ft]). Then (|70|) takes the form 

50 / 1 1 \ 

dk f{E, k) . + MJ? ^ , (73) 

where the separate contributions of electrons and positrons are evident. If the 
poles are second-order we have to take care of the contribution coming from the 
derivative of the ^-function integrated through the C± contours. This generates 
the terms 

/OO 1 POO 1 

^ dk h e > k K^) + i S{E ~" ) + L„ dk HE - k K-^ + i S{E - = 

K§)2/ (£ ^7^^ /(£ ^ ) ' (74) 



This is a T independent term. When inserted in ( |70D it cancel the result of the 
integral over C . 
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